Advances in point-contact spectroscopy: two-band superconductor MgB 2 (A review) 
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Analysis of the point-contact spectroscopy (PCS) data on the new dramatic high-T c supercon- 
ductor MgB2 reveals quite different behavior of two disconnected a and n electronic bands, deriving 
from their anisotropy, different dimensionality, and electron-phonon interaction. PCS allows direct 
registration of both the superconducting gaps and electron-phonon-interaction spectral function 
of the two-dimensional a and three-dimensional tt band, establishing correlation between the gap 
value and intensity of the high-T c driving force - the Ei g boron vibration mode. PCS data on some 
nonsuperconducting transition-metal diborides are surveyed for comparison. 
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MgP>2 was discovered to be superconducting only a 
couple of years ago pj, and despite that many of its 
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characteristics have now been investigated and a con- 
sensus exists about its outstanding properties. First of 
all, this refers to its high T c 40 K) which is a record- 
breaking value among the s-p metals and alloys. It ap- 
pears that this material is a rare example of multi-band 
(at least two) electronic structure, which are weakly con- 
nected with each other. These bands lead to very un- 
common properties. For example, T c is almost indepen- 
dent of elastic scattering, unlike for other two-band su- 
perconductors 0- The maximal upper critical magnetic 
field can be made much higher than that for a one-band 
dirty superconductor The properties of MgB 2 have 
been comprehensively calculated by the modern theoreti- 
cal methods, which lead to a basic understanding of their 
behavior in various experiments. 



Crystal structure 

Magnesium diboride, like other diborides MeB2 
(Me=Al, Zr, Ta, Nb, Ti, V etc.), crystalizes in a hexag- 
onal structure, where honeycomb layers of boron are in- 
tercalated with hexagonal layers of magnesium located 
above and below the centers of boron hexagons (Fig. 
The bonding between boron atoms is much stronger than 
that between magnesium, and therefore the disordering 
in the magnesium layers appears to be much easier than 
in the boron layers. This difference in bonding between 
boron and magnesium atoms hinders the fabrication of 
MgB 2 single crystals of appreciable size. 



Electron band structure 

The electron band structure of MgB 2 has been calcu- 
lated using different ab initio methods yielding basically 
the same result 0, HUES 13 • The E(k) curves are shown 
in Fig.EJ 

The dispersion relations are shown for boron p- 
character orbitals, which play a major role in trans- 
port and thermodynamic properties. The radii of the 
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FIG. 2: Band structure of MgB2 with the B p-character. The 
radii of the hollow (filled) circles are proportional to the ir (a) 
character and zero-line marks the Fermi energy. After Mazin 
et aZ.0. 

hollow circles are proportional to the 7r-band character, 
which is made from p z boron orbitals, while those of 
the filled circles are proportional to the cr-band charac- 
ter, made from p xy orbitals. The most important is a 
quasi two-dimensional dispersion relation along the TA 
(A) direction with a small Fermi energy «0.6 eV, and 
accordingly, with a moderate Fermi velocity. The corre- 
sponding sheets of the Fermi energy form the cylindri- 
cal surfaces along the TA direction seen in Fig.[S] below. 
The corresponding electron transport is very anisotropic 
(p c / Pab — 3.5 0|) with the plasma frequency for the a 
band along the c (or z) axis being much smaller than 
that in the ab (xy) direction . The hole branch along 
TA experiences a huge interaction with the phonon E 2g 



mode for carriers moving along the ab plane (see below) , 
although its manifestation is screened effectively by the 
much faster hole mobility in the 7r-band Q. 

In a dirty material, with prevailing disorder in the 
magnesium planes, the 7r-band conductivity is blocked 
by defects, and the a band takes over, implying greater 
electron phonon interaction (EPI) than in the clean ma- 
terial. This constitutes a plausible explanation for the 
violation of the Matthiesen rule, which manifests itself 
in an increase of the residual resistivity together with an 
increase of the temperature coefficient at a high temper- 
ature @- 

At the same time, the critical temperature T c does not 
decrease substantially in dirty materials Q , since the su- 
perconductivity is induced by EPI in the cr-band, whose 
crystal order is much more robust. 

This consideration is very important in understanding 
the point-contact data, since the disorder at the surface of 
the native sample depends on the position of the contact 
spot, and because of the uncontrolled introduction of the 
further disorder while fabricating the contact. 



Critical magnetic field 

In a clean material the layered crystal structure dic- 
tates strong anisotropy of the upper critical magnetic 
fields S°| 3> B c c2 . Their ratio at low temperatures 
reaches about 6 while B^ 2 is as low as 2-3 T jj^- If 
the field aligned is not precisely along the ab plane, the 
B c2 value is strongly decreased. 

On the other hand, for a dirty material the anisotropy 
is decreased (to a ratio of about 1.6—2), but both the 
magnitudes of B^ 2 and B c c2 are strongly increased. For 
strongly disordered sample, it may be as high as 40 T 
0! It is interesting that this high value is achieved at 
low temperature, where the disordered 7r band is fully 
superconducting. 

Hence, we may expect that the value of critical mag- 
netic field at low temperatures is the smaller the cleaner 
is the part of the MgB2 volume near the contact, pro- 
vided its T c ~ yowK, This observation is important in 
the classification of contacts with respect to their purity. 



Phonons and Electron-Phonon Interaction 

The phonon density of states (PDOS) is depicted in 
Fig. The upper panel shows the measured PDOS at 
T = 8K, while the lower ones shows the calculated DOS 
with the partial contribution from boron atoms moving 
in the a6-plane and out of it. One can see the peak for 
boron atoms moving in the ab plane at ~ 75 meV, which 
plays a very important role in the electron-phonon in- 
teraction, as is shown in FigJU measured by inelastic 
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FIG. 3: Upper panel: Phonon density of states in MgB2 
determined experimentally by neutron scattering. Bottom 
panel: calculated curve (solid line) with decomposition on 
boron atoms vibrating out of ab plane (dotted curve) and 
parallel to it (dashed curve). After Osborn et aZ.|l3j|. 



X-ray scattering This mode gives a weakly disper- 
sion branch between 60 and 70 meV in the TA direction 
with E2 g symmetry at the V point. The linewidth of 
this mode is about 20-^28 meV along the TA direction, 
while along the TM direction it is below the experimen- 
tal resolution. The same phonon peak is active in Raman 
scattering fl5t Ht| . It is located at the same energy 
with the same linewidth. This points to the very strong 
EPI for this particular lattice vibration mode. The same 
result follows from theoretical considerations. 

Figure [3] shows the distribution of the superconducting; 
energy gap on the Fermi surface of MgB2 |lj| ■ The max- 
imum gap value is calculated along the TA direction due 
to the very strong EPI. Just in this direction is located 
2D a band (cylinders along the TA direction). The 3D 
7r band has a much smaller EPI, and, correspondingly, a 
smaller energy gap. The EPI parameter A can be decom- 
posed between different pieces of the Fermi surface. It is 
shown 0| that the value of A on the a band amounts to 
2 — 3. Moreover, X a can be decomposed between different 
phonon modes, and it appears that only the E2 ff phonon 




FIG. 4: Dispersion curves of phonons in MgB2 and the width 
of phonon lines determined by inelastic X-ray scattering (sym- 
bols) together with calculations (solid lines). After Slmkla 
et aZ.pf. 

mode along the TA direction plays a major role with a 
partial \ a value of about ~ 25 |20j, though concentrated 
in a very restricted phase space. 

Mechanism for high T c in MgB2 

The commonly accepted mechanism for high T c in 
MgB2 is connected with the strong interaction between 
charge carriers and phonons in the Ei g mode. This mode 
is due to antiparallel vibration of atoms in the boron 
planes. The key issue is that along the TA direction the 
electron band structure is such that the Fermi energy of 
the hole carriers is only 0.5-f-0.6 eV, which shrinks even 
more when the borons deviate from the equillibrium posi- 
tions. Together with the 2D structure of the correspond- 
ing sheet of the Fermi surface, this leads to a constant 
density of states at the Fermi energy and, correspond- 
ingly, to very large EPI with partial \„ (the EPI pa- 
rameter in the a band) of about ~ 25 [2rj. Cappelluti 
et al. 21] point out that the small Fermi velocity for 
charge carriers along the TA direction leads to a large 
nonadiabatic correction to T c (about twice as much com- 
pared with the adiabatic Migdal-Eliashberg treatment). 
Although this interaction is a driving force to high T c 
in this compound, it does not lead to crystal structure 
instability, since it occupies only a small volume in the 
phase space. 

The role of the ir band is not completely clear. On the 
one hand, the tt and a bands are very weakly connected, 
and for some crude models they can be thought as be- 
ing completely disconnected. On the another hand, the 
energy gap of the 7r band goes to zero at the same T c 
as in the bulk, and correspondingly 2A 7r (0)/fcT c = 1.4, 
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FIG. 5: Superconducting energy gap distribution over the 
Fermi surface (FS) of MgB2. The gap value around 7 meV 
corresponds to cylinder like sheets of the FS centered at V 
points, while the small gap value around 2 meV corresponds 
to the tubular FS network. After Choi et al.UM. 



which is much less than the value predicted by the weak 
coupling BCS theory. One can think of the ir band as 
having intrinsically much lower T c 10 K than the bulk 
[22| . and at higher temperatures its superconductivity is 
induced by a proximity effect in k-space from a band 
\2?\ . This proximity effect is very peculiar. On the one 
hand, this proximity is induced by the interband scatter- 
ing between the it and a sheets of the Fermi surface. On 
the other hand, the charge carriers connected with the 
7r band are mainly located along the magnesium planes, 
which can be considered as a proximity effect in coordi- 
nate space for alternating layers of S — N — S structure, 
although on a microscopic scale. Moreover, many of the 
unusual properties of MgB2 may be modelled by an al- 
ternating S — N — S layer structure, the limiting case 
to the crystal structure of MgB2. In other words, MgB2 
presents a crossover between two-band superconductivity 
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FIG. 6: Scanning electron microscopy image of MgB2 films. 
After Kang et aZ^2f|. 



and a simple proximity effect structure. 



SAMPLES 

We have two kind of samples supplied for us by our 
colleagues from the Far-East. [S^ 

The first is a thin film with a thickness of about several 
hundred of nanometers (Fig.[5J) [24|. Similar films have 
been investigated by several other groups with different 
methods. These films are oriented with their c-axis per- 
pendicular to the substrate. The residual resistance is 
about several tens of /ifl cm with a residual resistance 
ratio (RRR) ~ 2.2. This means that on average the films 
have a disorder between crystallites. 

It does not exclude the possibility that on some spots 
the films contain clean enough small single crystals on 
which we occasionally may fabricate a point contact: see 
Fig. [fjj Normally, we make a contact by touching the film 
surface by noble metal counter electrode (Cu, Au, Ag) 
in the direction perpendicular to the substrate. Thus, 
nominally the preferential current direction in the point 
contact is along the c axis. Nevertheless, since the sur- 
face of the films contains terraces with small crystallites, 
point contact to the ab plane of these crystallites is also 
possible. Sometimes, in order to increase the probability 
of making the contact along the ab plane, we broke the 
substrate with the film and made contact to the side face 
of the sample. 

The second type of sample is single crystal [jfij which 
also was measured by other groups [nj, l27| . Crystals 
are plate-like (flakes) and have sub-millimeter size (see 
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FIG. 7: Scanning electron microscopy image of MgB2 single 
crystals. After Lee el al. |26| . 



Fig. 01 . They were glued by silver epoxy to the sample 
holder by one of their side faces. The noble metal counter 
electrode was gently touched in liquid helium by another 
(the opposite) side face of the crystal. In this way we try 
to preferentially make a contact along the ab plane. On 
average, in the bulk, the single crystals are cleaner than 
the films, but one should be cautious, since the properties 
of the crystal surface differ from the properties of the 
bulk, and fabrication of a point contact may introduce 
further uncontrolled defects into the contact area. 

Thus, a priori one cannot define the structure and com- 
position of the contacts obtained. Nevertheless, much of 
that information can be ascertained by measuring vari- 
ous characteristics of a contact. Among those the most 
important is the Andreev-reflection-non-linearities of the 
I — V curves in the superconducting energy-gap range. 
The magnetic field and temperature dependences of the 
superconducting non-linearities supply us with additional 
information. And finally, much can be extracted from 
the I — V nonlincarities in the normal states (the so- 
called point-contact spectra). The more information we 
can collect about the electrical conductivity for different 
conditions of the particular contact, the more detailed 
and defined picture of it emerges. It is not an easy task, 
since a contact has limited lifetime, due to electrical and 
mechanical shocks. 

Let us make a rough estimate of the distance scales 
involved in the problem. The crystallite size of films is of 
the order of 100 nm (see H^). The contact size d in the 
ballistic regime equals d ~ \J pl/R (the Sharvin formula) . 
Taking pi ^ 7 10~ 7 f7 cmx7 10" 6 cm= 5xl0- 12 ficm 2 (lfij . 
we obtain d ~ 7 nm both along the ab and c directions 
for typical resistance of 10 £1. If we suppose that a grain 
is dirty (with very short mean free path), then we apply 



the Maxwell formula d ~ p/R with the results for d about 
0.7 nm and 2.6 nm for ab and c directions, respectively, 
taking p for corresponding directions from the same ref- 
erence ^c|- Thus, the contact size can be of the order 
or smaller than the electronic mean free path (l a t = 70 
nm and l c — 18 nm, according to |lfj|). which means that 
we are working, apparently, in the spectroscopic regime, 
probing only a single grain. 

Rowell |28| . analyzing a large amount of experimental 
data for the resistivity and its temperature dependence, 
came to the conclusion that for highly resistive samples 
only a small part of the effective cross section should be 
taken into account. The reason is that the grains in MgB2 
are to great extent disconnected by oxides of magnesium 
and boron. For point-contact spectroscopy previous anal- 
ysis leads us to the conclusion that the contact resistance 
is frequently measured only for a single grain or for sev- 
eral grains, with their intergrain boundaries facing the 
contact interface. This is due to the current spreading 
on a scale of the order of the contact size d near the 
constriction. 



THEORETICAL BACKGROUND OF PCS 
Non-linearity of I — V characteristic 

The non-linearities of the I — V characteristic of a 
metallic contact, when one of the electrodes is in the 
superconducting state, can be written as |22, 30] 

I (V) c ^ - 5I^ h (V) + I exc (V) (1) 

Here Rq is the contact resistance at zero bias in the nor- 
mal state. 8Ip h (V) is the backscattering inelastic current 
which depends on the electron mean free path (mfp) I. 
For ballistic contact this term is equal, in order of mag- 
nitude, to 



si&(y) ~ -i(v) (2) 

Hn 

where U n is the inelastic electron mfp, and d is the char- 
acteristic contact diameter. If the electron flow through 
the contact is diffusive (l e i <C d, l e i being an elastic mfp) 
but still spectroscopic, since \fU n l e i ^> d, then the ex- 
pression J5J) should be multiplied by l e i jd. This decreases 
the characteristic size, for which the inelastic scattering 
is important, from d to l e \ (d — ► l e i), and for short l e i 
makes the inelastic current very small. We notice that 
the inelastic backscattering current SIp h (V) in the super- 
conducting state is approximately equal to the same term 
in the normal state. Its second derivative turns out to 
be directly proportional to the EPI function a 2 (oj) F(oj) 
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dV 2 3 hv F 

where a describes the strength of the electron interaction 
with one or another phonon branch, and F(lu) stands 
for the phonon density of states. In point-contact (PC) 
spectra the EPI spectral function a 2 (uj) F(cu) is modi- 
fied by the transport factor, which strongly increases the 
backscattering processes contribution. 

In the superconducting state the excess current 
lexc 0; which is due to the Andreev reflection of electron 
quasiparticles from the N-S boundary in a N-c-S point 
contact (c stands for " constriction" ) , can be written as 



/exc (V) = I° exc + 5I exc (V) (4) 

where I® xc ~ A/Rq ss const for eV > A (A being the 
superconducting energy gap). 

The nonlinear term in the excess current Q can be 
decomposed in its turn in two parts, which depend in 
different ways on the elastic scattering of electron quasi- 
particles: 
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{V) = Slf xc {V) + 5I™ c {V) (5) 

where 51* e (V) is of the order of (A/eV)I° xc , and 
SI ixc( v ) ~ Wkn)Iexc- Notice that the latter be- 
haves very similar to the inelastic backscattering current 
SIp h (V), namely, it disappears if l e i — > 0, while the first 
term in the right hand side of expression JSJ does not de- 
pend on l e i in the first approximation. This enables one 
to distinguish the elastic term from the inelastic. Finally, 
all excess current terms disappear when the superconduc- 
tivity is destroyed, while 8I^ h {V) remains very similar in 
both the superconducting and normal states. 

The expression for the elastic term in the excess 
current was calculated for ballistic N-c-S contacts by 
Omelyanchuk, Kulik and Beloborod'ko |33J. Its first 
derivative equals (T = 0) : 



dlt 



dV 



ballistic 



NcS 



1 



A(eV) 



zV + yj (eV) 2 - A 2 (eV) 



(6) 



For the diffusive limit (k <C d) , Beloborod'ko and Ku- 
lik derived the current-voltage characteristic (see Eq.(21) 
in Ref.[34j]), which for the first derivative at T = gives 
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dV 
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eV + A(eV) 
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eV 



v /(eV) 2 -A 2 (eV) 



A(eV) 



V(eV) 2 -A2(eV) 



(7) 



For the sake of comparison, the similar expression of 
the nonlinear term in NIS tunnel junctions (I stands for 
"insulator"), due to the self-energy superconducting en- 
ergy gap effect, has the form [3fij ]: 



dl\ 1 m 

— = — 

dv J nis R o 



eV 



{eVf - A 2 [eV) 



(8) 



Equations , Q , and JSJ are identical in their structure 
and take into account the same effect, viz., the renormal- 
ization of the energy spectrum of a superconductor in the 
vicinity of characteristic phonon energies. 

From the expressions (0), (@J, and (JSJ it becomes 
clear that only on the relatively clean spots one can ob- 
serve the inelastic backscattering current SIp h (V) pro- 
vided that the excess current term 8I" XC (V) is negligible. 
The latter can be canceled by suppression of supercon- 
ductivity either with magnetic field or temperature. On 
the contrary, in the superconducting state, for dirty con- 
tacts, all the inelastic terms are very small, and the main 
nonlinearity is provided by the A(ey)-dependence of the 
excess current J7J). 



Two-band anisotropy 

Brinkman et al. have shown that in the clean case 
for an NIS MgB2 junction, the normalized conductance 
is given by 

„ m (ft) nis H) 2 aAV) + H) 2 °AV) (q) 

\ ) Tdr\ ; — n — ; — n l y ) 



\dVJNIN 



where 



•n-(cr) 



is the plasma frequency for the 7r (a) band 



and ^ (V) is the normalized conductivity of the ir(a) 
band separately. The calculated tunneling conductance 
in the afe-plane and along the c-axis are |l 11 

a ah (V) = 0.67a, (V) + 0.33a, (V) (10) 
a c (V) = 0.99OV (V) + 0.01(T a (V) (11) 

Hence, even along the aft-plane the contribution of the 
a band is less than that of the tt band, to say nothing 
about the direction along the c axis, where it is negligible 
small. The calculation predicts that if the " tunneling 
cone" is about several degrees from the precise ah plane, 
then the two superconducting gaps should be visible in 
the tunneling characteristics. In other directions only a 
single gap, corresponding to the 7r band, is visible. We 
will see below that this prediction is fulfilled in a point- 
contact experiment, as well. 

Things are even worse when one tries to measure the 
anisotropic Eliashberg function by means of supercon- 
ducting tunneling. The single-band numerical inversion 
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FIG. 8: Typical shapes of dV/dl (experimental dots) for 4 
contacts between MgB2 thin film and Ag with the correspond- 
ing BTK fitting (lines) |4C||. &l{8) stand for large (small) 
superconducting energy gap. After Naidyuk et al. 40]. 



program |3fil Is?)] gives an uncertain result, as was shown 
in Ref. jaJT 

Point-contact spectroscopy in the normal state can 
help in this deadlock situation. It is known that the in- 
elastic backscattering current is based on the same mech- 
anism as an ordinary homogeneous resistance, provided 
that the maximum energy of the charge carriers is con- 
trolled by an applied voltage. The electrical conductivity 
of MgB2 can be considered as a parallel connection of two 
channels, corresponding to the ir and a bands .2]. The 
conductivity of the 7r band can be blocked by disorder of 
the Mg-atoms. This situation is already obtained in ex- 
periment, when the temperature coefficient of resistivity 
increases simultaneously with an increase of the residual 
resistivity, which leads to violation of Matthiessen's rule 
(see Fig. 3 in 0). In this case we obtain the direct access 
to the cx-band conductivity, and the measurements of the 
PC spectra of the EPI for the a band is explicitly possible 
in the normal state. Below we will see that this unique 
situation happens in single crystals along ab plane. 



EXPERIMENTAL RESULTS 
Superconducting energy gaps 

c-axis oriented thin films 

Our measurements of the superconducting energy gap 
by means of Andreev reflection from about a hundred 
NcS junctions yield two kinds of dV/dl-cmve, shown in 
Fig-El The first one clearly shows two sets of energy 
gap minima located, as shown in distribution graph of 
Fig.El(upper panel), at 2.4± 0.1 and 7.1±0.4 meV. These 
curves can be nicely fitted by BTK |3^| theory (with 
small r parameter) for two conducting channels with an 
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FIG. 9: Superconducting energy gap distribution of about 100 
different junctions prepared on a MgB2 c-axis oriented film. 
On the lower panel the theoretical distribution is shown. After 
Naidyuk et aZ.^}. 



adjusted gap weighting factor [4fJ. The second kind is 
better fitted with a single gap provided an increased de- 
pairing parameter T (Fig. [5] (middle panel)). Certainly, 
the division of the gap structure into two kinds mentioned 
is conventional, and depends upon the circumstance that 
the larger energy gap is explicitly seen. These two kinds 
of gap structure comprise about equal parts of the to- 
tal number of junctions. Usually the contribution of the 
large gap in the double-gap spectra is an order of mag- 
nitude lower than that of the small one, which is in line 
with the small contribution of the a band into the con- 
ductivity along the c axis (see Ea,ITT|l. 

It is important to note that the critical temperature of 
the material around the contact is not more than a few K 
below T c in the bulk material. This is determined by the 
extrapolating the temperature dependence of PC spec- 
tra up to the normal state. Such an insensitivity of T c on 
the elastic scattering rate is explained in Ref. [2| . Never- 
theless, we stress that the gap structure (either double- 
or single-gap feature, and the position of the single-gap 
minimum on dV/dl) depends very much on random vari- 
ation of the scattering in the contact region. Moreover, 
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FIG. 10: Negative slope of dV/dl at large biases for a 36 
fi contact between MgB2 single crystal and Ag showing the 
magnetic-field gap-structure evolution at 4.2 K. 



FIG. 11: Large gap structure evolution for single crystal 
MgB2-Au 87 Q junction in magnetic field at 4.2 K. The curves 
are shifted vertically for clarity. 



since the main part of the junction conductivity is due to 
the charge carriers of the n band, even the background 
conductance quite often follows the " semiconductive" be- 
havior, namely, the slope of the dV/dl curve at large bi- 
ases is negative (Fig. I10[l . That means, that the carriers 
in the tt band are close to localization |4lJ . 

In the lower panel of Fig.[5] the theoretical prediction 
of the energy gap distribution |18( is shown. One can 
see that the theoretical positions of the distribution max- 
ima coincide approximately with the experimental values. 
Only the low-lying maximum is not seen in the experi- 
ment. It should be noted that according to Mazin et al. 
|42| variation of the superconducting gaps inside the a 
and 7T bands can hardly be observed in real samples. 

The distribution of the different gaps over the Fermi 
surface is shown in Fig. |3J One can immediately see that 
for a c-axis oriented him the main structure should have 
a smaller gap, which is approximately isotropic. Only if 
the contact touches the side face of a single crystallite 
(Fig. EJ), is the larger gap visible, since it corresponds 
to the cylindrical parts of the Fermi surface with Fermi 
velocity parallel to the ab plane. 



Single crystals 

The same variety of energy gap structure is observed 
for single crystals as well, but with some peculiarity due 
to preferential orientation along the ab plane. The most 
amazing of them is the observation of dV/dl-g&p struc- 
ture in Fig. ^2 with visually only the larger gap present. 
This gap persist in a magnetic field of a few tesla un- 
like the smaller gap, which according to (Islliij vanishes 
above 1 T. Spectra of that kind were not observed in thin 
films. This means that the conductivity is governed only 
by the a band. This may be caused by the circumstance 
that the tt band is blocked completely by Mg disorder or 



by oxidation of Mg atoms on a6-side surface of the crys- 
tal. At the same time, in a single crystal there is much 
less scattering in the boron planes, due to the robustness 
of the B-B bonds. We will see below that just this case 
enables us to observe directly the most important Ei g 
phonon mode in the electron-phonon interaction within 
the a band. 

In single crystals the negative slope in dV/dl curve at 
large biases is observed quite often, which confirms that 
the disorder in the tt band leads to quasi-localization of 
charge carriers. An example of this is already shown in 

Fig. uni 

FigureslT^landlTSl display a series of magnetic-field and 
temperature dependence of dV/dl curves with their BTK 
fit. Here the two gaps are clearly visible, corresponding 
to the theoretical prediction in the ab direction Eq. 1)11(1 . 
The temperature dependence of both gaps follows the 
BCS prediction (see Fig. ll4|) . For temperatures above 
25 K their behavior is unknown because this particular 
contact did not survive the measurements likely due to 
thermal expansion of the sample holder. 

Figure ED displays the magnetic field dependences of 
large and small gaps. Surprisingly, the small gap value is 
not depressed by a field of about 1 T, and the estimated 
critical field of about 6T is much higher as stated in 
0, 0| , although the intensity of the small-gap minima 
is suppressed rapidly by a field of about 1 T. Correspond- 
ingly, the small-gap contribution w 56] to the dV/dl 
spectra is decreased by magnetic field significantly, from 
0.92 to 0.16 (see Fig. ll5|) . while w versus temperature 
even increases slightly from 0.92 at 4.3 K to 0.96 at 24 K 
(not shown). 

The area under the energy-gap minima in dV/dI(V) 
is approximately proportional to the excess current I exc 
(see Eq.(@J) at eV ^> A (or roughly to the superfluid den- 
sity). The excess current depends on the magnetic field 
with a positive overall curvature fFig. H6fl . I exc (B) de- 
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FIG. 12: Magnetic field dependences of dV/dl curves (solid 
lines) for a single crystal MgB2-Cu 2.2S7 junction along the 
ab plane with their BTK fittings (thin lines). Two separate 
sets of gap minima are clearly seen at low fields. 
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FIG. 13: Temperature dependences of dV/dl curves (solid 
lines) for the same junction as in Fig. ll2l with their BTK fit- 
tings (thin lines). 




FIG. 15: Magnetic field dependences of the large and small su- 
perconducting energy gaps (solid triangles) obtained by BTK 
fitting from Fig. 1121 Open triangles show the F value for large 
and small gap, respectively. The circles demonstrate the de- 
pression of the small-gap contribution to the dV/dl spectra 
by magnetic field. The lines connect the symbols for clarity. 




FIG. 16: I exc (B) (squares) for a MgB2-Cu junction from 
Fig. 1121 The dashed lines show the different behavior of 
Iexc{B) at low and high fields. 



6 - 




1 1 ' 1 ' 1 ' 1 ' 1 ' ' 1 1 

5 10 15 20 25 30 

T(K) 

FIG. 14: Temperature dependences of large and small super- 
conducting energy gaps obtained by BTK fitting from Fig. 1131 
The solid lines represent BCS-like behavior. 



creases abruptly at first and then more slowly above 1 T. 
This corresponds to a drastic depressing of the dV/dI(V) 
small-gap-minima intensity by a magnetic field of about 
1 T and to robustness of the residual superconducting 
structure against further increase of magnetic field. This 
is a quite different dependence from what is expected for 
Iexc, which is in general proportional to the gap value 

Q. 

In contrast, I e xc(T) has mostly negative curvature and 
shape similar to the BCS dependence. Often a positive 
curvature appears above 25 K (see, e.g. Fig. lT%)l . 

This kind of anomaly can be due to the two-band 
nature of superconductivity in MgB2, since the mag- 
netic field (temperature) suppresses the superconductiv- 
ity more quickly in the ir band and then, at higher field 
(temperature), in the a band. The same consideration 
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FIG. 17: Temperature dependence of the penetration depth 
in the model of two independent BCS superconducting bands 
(dashed and dotted line) with different superconducting gaps. 
The resulting penetration depth (solid line) clearly shows a 
non-BCS temperature behavior. The low-temperature behav- 
ior will be dominated by the band with the smaller supercon- 
ducting gap. After Golubov et al. |46j . 




T (K) 

FIG. 18: Temperature dependence of a single superconduct- 
ing energy gap (squares) obtained by BTK fitting of dV/dl 
curves from inset. The solid lines represent BCS-like behav- 
ior. The triangles show the dependence of the excess current. 
Inset: dV/dl curves for a MgB2-Cu 8il contact at different 
temperatures. 



is valid for 1/Xl, which is roughly proportional to the 
" charge density of superfiuid condensate" . In the case of 
zero interband scattering, the simple model |46| predicts 
the temperature dependence shown in Fig. 1171 for a and 
7r parallel channels, which will yield a smooth curve with 
general positive curvature, taking into account the small 
interband scattering occurring in reality. 

If the 7r-band conductivity is blocked by a short mean 
free path, then the curvature of I exc (T), being propor- 
tional to A CT (T), should be negative, which supplies us 
with additional confirmation of single band conductivity 
along the a band. Thus, measuring the magnetic field 
and temperature dependences of I exc can elucidate the 
contact structure. 

Figure IT81 displays the temperature dependence of the 
gap for the dV/dl curves with a single gap structure, 
which vanishes around 25 K. A magnetic field of 1 T sup- 
presses the gap minima intensity by factor of two, but 
the minima are clearly seen even at 4T (not shown), the 
maximal field in this experimental trial. This excludes an 
origin of these gap minima due to small gap. According 
to the calculation in |ll| a large amount of impurity scat- 
tering will cause the gaps to converge to A ~4.1 meV and 
T c to 25.4 K. Therefore these single-gap spectra reflect 
a strong interband scattering due to impurities, which 
likely causes a "semiconducting- like" behavior of dV/dl 
above T c (see Fig.^J inset). I exc (T) behaves nearly as 
A(T) except in the region T>25K, where I exc is still 
nonzero because of a residual shallow zero-bias minimum 
in dV/dl above 25 K. 
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FIG. 19: Raw PC EPI spectrum for a ZrB2 5.5 £1 point contact 
at 4.2 K. After Naidyuk et al$\. 



Phonon structure in the I — V characteristics 

PC EPI spectra of nonsuperconducting diborides 

We have studied the PC EPI spectra d 2 V/dI 2 cx 
—d 2 I/dV 2 (see also Eq. Q) of non-superconducting di- 
borides MeB 2 (Me=Zr, Nb, Ta) 47]. The cleanest sam- 
ple we have is a ZrB 2 single crystal, and its PC EPI 
spectrum is shown in Fig. 1191 One recognizes a classical 
PC EPI spectrum from which one can estimate the posi- 
tion of 3 main phonon peaks and obtain the lower limit 
of EPI parameter A pc 47] . 

Essentially similar spectra were observed for other di- 
borides, taking into account their purity and increased 
EPI, which leads to a transition from the spectroscopic to 



11 





\ ZrB 2 ■ 




% 







ij 
y'.i 


■j 











100 

Energy (meV) 



FIG. 20: Comparison of high-resolution electron-energy- 
loss spectroscopy measurements of surface phonon dispersion 
(bottom panels, symbols) [ifl^ with the PC spectra for ZrB2 
and NbP>2 after subtraction of the rising background (upper 
panels) . 



a non-spectroscopic (thermal) regime of current flow [42fl • 
The positions of the low-energy peaks are proportional 
to the inverse square root of the masses of the d met- 
als ^?J' as expected. For these compounds the phonon 
density of states is measured by means of neutron scat- 
tering |48| and the surface phonon dispersion is derived 
by high-resolution electron-energy-loss spectroscopy . 
The positions of the phonon peaks or duj/dq=0 for the 
dispersion curves correspond to maxima of the PC spec- 
tra (Fig.EOI Fig.EH). 



PC EPI spectra of MgB-z in c-axis oriented films 

From the above consideration we had anticipated that 
one could easily measure the EPI spectral function of 
MgB2 in the normal state, provided that the supercon- 
ductivity is destroyed by magnetic field. Unfortunately, 
that was not the case. The stronger we suppress the su- 
perconductivity in MgB 2 , the less traces of phonon struc- 
ture remain in the / — V characteristic and its derivatives 
(Fig. l22[) 2j^. This is in odd in relation to the classical 
PCS, since the inelastic phonon spectrum should not de- 
pend on the state of electrodes in the first approximation 
(see section Theoretical background). 

Instead, most of the MgB 2 spectra in the supercon- 
ducting state show reproducible structure in the phonon 
energy range (Fig. l23|) which was not similar to the ex- 
pected phonon maxima superimposed on the rising back- 
ground. This structure disappears upon transition to the 
normal state. Quite interestingly is that the intensity of 
this structure increases with increase of the value of the 
small gap, which means that the gap in the it band and 
observed phonon structure is connected j^. Based on 
the theoretical consideration mentioned above, we con- 
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FIG. 21: Comparison of phonon DOS neutron measurements 
after Heid et al. [48| (symbols) with PC spectra for NbB2 and 
TaB2 after subtraction of the rising background (solid curves) . 
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FIG. 22: Phonon singularities in the PC spectra of a MgP>2 
thin-film - Ag junction as a function of magnetic field and 
temperature. T and B are shown beside each curve. After 
Yanson et al. |23| . 



elude that the disorder in the it band is so strong that 
it precludes observation of the inelastic current, and the 
phonon nonlinearities of the excess current (see Eq.©) 
play the main role, which does not depend on the scat- 
tering. 

Very rarely we observed signs that the measured char- 
acteristics indeed satisfy the conditions imposed on the 
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FIG. 23: Superconducting gap features (upper panel) and 
phonon structure (bottom panel) in the spectra of thin-film 
MgB2-Ag junctions with different resistances at T = 4.2K, 
5 = 0. After Yanson et al. 23]. 



inelastic PC spectra. One such example is shown in 
Figs.|2J[5Bl For this particular junction the supercon- 
ducting peculiarities are almost completely suppressed 
above 20 mV by moderate field (4T). What remains is 
a weak zero-bias minimum (~ 1%) from the rather high 
value of the gap (Fig.^SJl. The background in dV/dI(V) 
rises nearly quadratically up to a few % of R at large 
biases (~100 mV). This leads to a linear background in 
V% oc d 2 V/dI 2 (V) with phonon peaks superposed above 
the background both in negative and positive bias polar- 
ity (compare with Fig. l*H?|l . The structure observed above 
30 meV corresponds reasonably in shape to the phonon 
density of states (Fig.|25J). At the low voltages (below 
30 meV), most probable, the gap peculiarities still pre- 
vail over the d 2 V/dI 2 (V) structure. Thus, for this con- 
tact we assume to observe the inelastic PC spectrum for 
the ir band, which should be compared to the Eliashberg 
EPI function for the same band calculated in Ref. | 51| 
(Fig l26fl . Both the experimental spectrum and the 7r- 
band Eliashberg function do not show anomalously high 
intensity of Ei g phonon mode, since only the Eliashberg 
function for a band is the principal driving force for high 
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FIG. 24: dV/dl and Vi oc d 2 V/dI 2 curves for a thin film 
MgB2-Ag junction revealing the inelastic PC spectrum for 
the 7r band. After Bobrov et al. |50j . 

T c in MgB 2 . The same conclusion should be ascribed to 
the excess-current phonon structure, since it also corre- 
sponds to the 7r band. This band has much larger Fermi 
velocity and plasma frequency along the c-axis compared 
to the a band [ll|. 

Thus, in order to register the principal EPI with the 
Ei g phonon mode, we are faced with the necessity of 
measuring the PC spectra for only the a band. This 
can be done in a single crystal along the ab plane with 
blocked 7r— band conductivity. 

PC EPI spectra of MgE>2 in the ab direction 

The desired situation was described in Ref. for sin- 
gle crystal oriented in the ab plane. As was mentioned 
above, the nominal orientation of the contact axis to be 
parallel to ab plane is not enough to be sure that this 
situation occurs in reality Moreover, even if one estab- 
lishes the necessary orientation (i. e., contact axis parallel 
to ab plane) the spectra should reflect both bands with 
a prevalence of the undesired ir band, because due to 
spherical spreading of the current the orientational selec- 
tivity of metallic point contact is much worse than that 
for the plane tunnel junction, where it goes exponentially. 
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FIG. 25: Comparison of the PC EPI spectrum (upper panel) 
from Fig. 1241 (after subtraction of the linear background and 
zero-bias maxima below 25 meV) with the phonon DOS mea- 
sured by neutron scattering || ^| (bottom panel) . 
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FIG. 26: Calculated Eliashberg functions for the a and 7r 
bands (inset). After Golubov et al. |5lf . 



The large mixture of 7r-band contribution is clearly seen 
from the gap structure in Fig. 1271 Inside the wings at 
the biases corresponding to the large gap (supposed to 
belong to the cr-band gap) the deep minima located at 
the smaller gap (correspondingly to the 7r-band gap) are 
clearly seen (see bottom panel of Fig. l27|) . The EPI spec- 



FIG. 27: Vi oc d 2 V/dI 2 (for two bias voltage polarities) 
and dV/dl curves for a single crystal MgB2-Cu junction 
(Ro = 1.5 fi) along the ab plane. Here the conductivity along 
the 7r band prevails, as is shown by the pronounced small-gap 
structure for the zero-field dV/dl curve at 4.5 K. The a 2 F(to) 
curve is the theoretical prediction for the 7r-band Eliashberg 
function from Fig. 1261 (inset) smeared similarly to the experi- 
mental data. After Naidyuk et al. 1521 . 



trum of the same junction is shown in the upper panel. 
One can see that the nonlinearities of the /— V character- 
istic at phonon biases are very small, and a reproducible 
structure roughly corresponding to the Eliashberg EPI 
function of the tt band [3a. l5*H appears in the bias range 
20 -i- 60 mV. Above 60 mV the PC spectrum broadens 
sufficiently hidden higher-lying phonon maxima. Even in 
the normal state (T > T c ), where the excess current dis- 
appears, one can see the kink at w 20 -f- 30 meV, where 
the first peak of the phonon DOS and the first maximum 
of the Eliashberg EPI function of the tt band occurs. At 
eV« 90 4 100 meV the PC EPI spectrum of Fig. ^sat- 
urates just where the phonon DOS ends. At T > T c in- 
termediate phonon peaks are hardly seen, since the ther- 
mal resolution, which equals 5.44 fc^T, amounts about 
20 meV, and the regime of current flow is far from ballis- 
tic, due to the high background observed. No prevalence 
of the E2 g phonon mode is observed, like a big maximum 
of EPI at w 60 4- 70meV or a kink at T > T c for these 
biases. 
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A quite different spectrum is shown in Fig.[5Sl which 
is our key result. Consider first the dV/dI{V) charac- 
teristics (see bottom panel). The energy gap structure 
shows the gap minima corresponding to the large gap 
(cr-band gap). The increase of dV/dI(V) at larger biases 
is noticeably larger than in the previous case fFig. 12711 . 
One can notice that the relatively small magnetic field 
(~1 T) does not decrease the intensity of gap structure 
substantially, unlike those for Fig. and even less than 
for Fig.EZI According to jH Ei] a field of about 1 T 
should depress the small gap intensity completely. All 
these facts evidence that we obtain a contact, in which 
only the cr-band channel in conductivity is operated. 

Let us turn to the PC EPI spectra d 2 V/dI 2 (V), which 
are connected via the following expression to the second 
harmonic signal V2 recorded in experiment: 



0,6 



1 d 2 V 
R 2 dl 2 



2^— 

V V 2 



Here R — dV/ dl and V\ is the rms value of the modula- 
tion voltage for the standard technique in the tunneling 
and point-contact spectroscopy. 

The PC EPI spectra for this contact are shown in 
Fig. 1281 (upper panel) for the highest field attainable in 
our experiments One can see that 8-^9 T is still 

not enough to destroy completely the superconductivity 
in the energy-gap low bias range (0-^30 meV), which can 
be taken as characteristic for a strongly superconducting 
a band. On the other hand, at larger biases no influ- 
ence of field is noted, which evidences that this part of 
I — V characteristic does not contains superconducting 
peculiarities, likely due to the high current density in the 
contact. Except for a small asymmetry, the spectrum 
is reproduced for both polarities. Before saturation at 
biases > 100 meV, where the phonon DOS ends, a well- 
resolved wide bump occurs, which is located at about 
60 meV. Further on, we will concentrate on this. 

First, we rescale it to the spectrum in R^dR/dV 
units, in order to compare with the theoretical estima- 
tion. We will show that the bump is of spectroscopic 
origin, i. e. is the regime of current flow through the con- 
tact is not thermal, although the background at large 
biases (V > 100 meV) is high. To do so, we compare 
this bump with a PC spectrum in the thermal regime for 
a model EPI function, which consists of a Lorentzian at 
60meV with small (2 meV) width. Calculated accord- 
ing to Kulik the thermal PC EPI spectrum is much 
broader, shown in Fig. 1291 as a dashed line. Any further 
increase of the width of the model spectra will broaden 
the curve obtained. Comparing the experimental and 
model spectra enable us to conclude that in spite of the 
large width, the maximum of the experimental spectra 
still corresponds to the spectroscopic regime. The high- 
temperature (T > T c ) spectrum in Fig.[2Hl shows the 
smeared kink at about 60meV, unlike that of Fig.1271 
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FIG. 28: Vi oc d 2 V/dI 2 (for two bias voltage polarities at 
9T) and dV/dl curves for a single-crystal MgB2-Cu junc- 
tion (Ro = 7.2 f2) along the ab plane. Here the conductivity 
along the a band prevails, as is shown by pronounced large- 
gap structure for the zero field dV/dl curve at 4.5 K. After 
Naidyuk et al.^. 
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FIG. 29: Comparison of the experimental spectrum of Fig. 1281 
with the thermal spectrum for a model spectral function in 
the form of a Lorentzian at 60meV with a width of 2meV 
(dashed line) and with the theoretical EPI spectra (bottom 
curves). After Naidyuk et al. 52]. 



15 




50 100 150 



1 1 I 1 


111,1111 
\ — 


--9T " 




1^1 . . . 


■■■ 1T 






— 0T 


1 


V tl0 mV 




i 



-100 100 

Voltage (mV) 



FIG. 30: Thermal limit for a band (as shows expressed large 
gap structure for zero field dV/dl curve at 4.5 K) in PC spec- 
trum of MgB2 single crystal along ab plane. After Naidyuk 
et al.^. 

Introducing greater disorder in the boron plane by a fab- 
rication procedure or by trying other spots on the side- 
face surface, the smeared thermal spectra were observed, 
coinciding in shape with the dashed curve in Fig.[3U] In 
this figure another junction is shown, where the energy 
gap structure also points to the u-band channel. Other 
junctions display the kink at about 30-^40 meV, like the 
high temperature spectrum in Fig. 1271 which together 
with their energy-gap structure can be ascribed to the 
thermal limit mainly in the ir band, despite the rather 
low bath temperature. 

A PC spectrum with broad maxima including also one 
at about 60 mV were observed in 0] on polycrystalline 
MgB2 samples driven to the normal state by applying 
moderate magnetic field and increasing of the tempera- 
ture. 

The large width of the EPI peak connected with the 
Ei„ phonon mode (Fig.l^SIl is not surprising. Shukla et 
al. [lj] measured the phonon dispersion curves along the 
TA and TM directions by means of inelastic X-ray scat- 
tering (see Fig.0J). The full width at half maximum for 
the E2 g mode along the TA direction amounts about 20- 
28 meV, which corresponds well to what we observe in 
the point-contact spectrum. If the phonon lifetime cor- 



responds to this (inverse) energy, then the phonon mean 
free path is about equal to the lattice constant and 
due to phonon reabsorption by accelerating electrons, we 
should anticipate a large background in the PC spectra 
as observed. If we compare the position of the bump 
(« 60 meV) with what is predicted for isotropic Eliash- 
berg EPI function 0] (see Fig. 1291) . then we, together 
with Shukla et al., should admit that the phonon-phonon 
anharmonicity is inessential for this mode, and its high 
width is due completely to the EPI. 

Now turn to the nonlinearity of the I — V curves due 
to electron-phonon interaction, which can be estimated 
from the dV/dl curves as about 10% for contact with 
the E 2s phonon modes in Fig.[2Hl This is comparable 
with the nonlinearity observed for nonsuperconducting 
diborides |47| with a small electron-phonon coupling con- 
stant. The reason for the relatively low nonlinearity of 
the I — V curves and low intensity of the principal E 2g 
phonon modes in the spectra for the MgB 2 contacts can 
be the fact that anomalous strong interaction is charac- 
teristic for restrictedgroup of phonons with sufficiently 
small wave vector |9(, whereas in point-contact spec- 
troscopy the large angle scattering is underlined. 

CONCLUSIONS 

We made an overview of the PCS investigations of c- 
axis oriented thin films and single crystals of MgB 2 . Our 
conclusions are as follows: 

1. There are two different superconducting gaps in 
MgB 2 , which are grouped at 2.4 and 7.0 meV. Roughly, in 
half of all point contacts studied for c-axis oriented films 
the two gap structure merges together due to strong elas- 
tic scattering remaining a single gap at about 3.5 meV. 

2. Anomalous temperature and especially magnetic 
field dependencies of excess current in point-contact junc- 
tions reflect the two-band structure of the superconduct- 
ing order parameter in MgB 2 . 

3. There are two mechanisms of revealing phonon 
structure in the PC spectra of MgB 2 : i) through the 
inelastic backscattering current, like for ordinary point- 
contact spectroscopy, and ii) through the energy depen- 
dence of the excess current, like in the similar tunneling 
spectroscopy of the electron-phonon interaction. They 
can be discriminated by destroying the superconductiv- 
ity with a magnetic filed and/or temperature, and by 
varying the electron mean free path. 

4. The prevailing appearance of the i? 2s boron mode, 
which mediates the creation of Cooper pairs, is seen in 
the PC spectra only along the a — b direction in accor- 
dance with the theory. The relatively small intensity of 
this mode in the PC spectra is likely due to their small 
wave vector and restricted phase volume. 

5. Related diborides (ZrB 2 , NbB 2 , and TaB 2 ) have 
d 2 V/dI 2 spectra proportional to the electron-phonon in- 
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teraction spectral function like that in common metals 
and a small EPI constant corresponding to their non- 
superconducting state. 
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Note added in proof 

After the paper was completed we have learned of the 
paper by Koshelev and Golubov where the mag- 
netic field dependence of A CT and A ff was presented. It 
turned out that the A a (B) and A T (i3) behavior is dif- 
ferent and is governed by diffusion constants depending 
on the coherence length. However, the critical field is 
the same both for A CT and A ff . This is in line with our 
observation given in Fig. 1151 Additionally, two experi- 
mental reports on the effect of magnetic field on both 
gaps in MgB2 by Gonnelli et al. (cond-mat/0308152) and 
Bugoslavsky et al. (cond-mat/0307540) appeared in the 
E-print archive. Bugoslavsky et al. reported that both 
order parameters persist to a common magnetic field. 
Gonnelli et al. corrected their previous claims and men- 
tioned that identification of the magnetic field at which 
the 7r-band features in dV/ dl visually disappear with the 
critical field for the tt band might not be correct. 
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